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Abstract. We address the question of the relations be- 
tween the black hole's mass, the accretion rate, the bolo- 
metric luminosity, the optical luminosity and the size of 
the Broad Line Region (BLR) in Active Galactic Nuclei, 
using recent observational data obtained from monitoring 
campaigns. We show that a standard accretion disc can- 
not account for the observed optical luminosity, unless it 
radiates at super-Eddington rates. This implies the exis- 
tence of another, dominant emission mechanism in the op- 
tical range, or a non standard disc (non stationary, ADAF 
and/or strong outflows). Narrow Line Seyfert 1 galaxies 
(NLSls) are most extreme in this context: they have larger 
bolometric to Eddington luminosity ratios than Broad Line 
Seyfert 1 (BLSls), and most likely a larger "non disc" 
component in the optical range. From realistic simulations 
of self-gravitating a-discs, we have systematically localized 
the gravitationally unstable disc and shown that, given 
uncertainties on both the model and observations, it co- 
incides quite well with the size of the BLR. We therefore 
suggest that the gravitationally unstable disc is the source 
which releases BLR clouds in the medium. However the 
influence of the ionization parameter is also required to 
explain the correlation found between the size of the BLR 
and the luminosity. In this picture the size of the BLR in 
NLSls (relative to the black hole size) is larger (and the 
emission line width smaller) than in BLSls simply because 
their Eddington ratio is larger. 
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an X-ray source, a photoionized Broad Line Region (BLR) 
and possibly a molecular torus at the scale of a few par- 
sees from the center (Antonucci & Miller 1985). Despite 
years of efforts, the physics of these components is far 
from being completely understood as well as the possible 
interaction between them. For instance, the origin of the 
observed correlation between the size of the BLR and the 
luminosity of the AGN, which, under some hypothesis, 
translates into a relation between the mass of the central 
black hole and the disc accretion rate, is still mysterious. 

In this paper, we try to understand the links between 
the accretion disc (at small and large radii), the optical, 
bolometric and Eddington luminosities, and finally the 
Broad Line Region (BLR). For this purpose, we use re- 
cent observational data which have led to the determina- 
tion of the size of the BLR, optical luminosity, and black 
hole (BH) mass in a few tens of objects spanning a large 
range of luminosity (Kaspi et al. 2000). In Section 2, we 
briefly recall the emission mechanisms in Active Galac- 
tic Nuclei (AGN), their link with the accretion disc, the 
main properties of the disc itself and those of the BLR. In 
Section 3 we discuss the relations between the observed 
optical luminosity, bolometric luminosity and BH mass. 
It is demonstrated that the standard disc emission cannot 
account for the optical luminosity, and that Narrow Line 
Seyfert 1 galaxies (NLSls) are the most extreme cases on 
this point of view. In Section 4, we model the outer disc 
and show that the occurrence of gravitational instability 
roughly coincides with the size of the BLR. Our conclu- 
sions are presented in the last section. 



1. Introduction 

There are several ubiquitous media in the central region 
of Active Galactic Nuclei (AGN): surrounding the super- 
massive black hole (BH), there are an accretion disc (AD), 
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2. General considerations 

2.1. Emission mechanisms in AGN: the inner disc 

Since its discovery in AGN spectra, the universal feature 
known as the "Big Blue Bump" (BBB) which represents 
the bulk of the bolometric luminosity, has been interpreted 
as thermal emission of an accretion disc (Shields 1978, 
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Fig. 1. Schematic view of the most central region of an AGN. The scaling depends slightly on the mass and on the 
accretion rate, and is appropriate for a 10 s Mq black hole (Rg — 1.5 X 10 13 cm) accreting at m ~ 0.1 in Eddington 
units. 



Malkan & Sargent 1982). If this disc is Keplerian and op- 
tically thick as generally assumed (the so-called "standard 
disc" model), then its emission reproduces roughly the 
BBB feature (see for a review Koratkar & Blaes 1999, and 
Collin 2000). Moreover, the idea that a non thermal power 
law continuum dominates the near infrared band and con- 
tributes, together with the disc, to the total emission in the 
optical range, has been abandoned after the publication 
of the generic spectrum of radio quiet quasars by Sanders 
et al. (1989). These authors have actually shown that the 
IR band is most likely due to hot dust emission peaking at 
a few microns. Since hot dust sublimates at temperatures 
~ 1700 K, dust cannot contribute to the optical emission. 
So, if one sustains the standard disc picture, one must 
admit that the optical continuum component is entirely 
emitted by the accretion disc. 

According to the standard model, the disc regions emit- 
ting in the optical band are located at a few 100 Rq, where 
Rq = GM/c 2 is the gravitational radius of the BH with 
mass M (see Fig. ^, middle panel) . Since the accretion lu- 
minosity varies as the optical luminosity should be 
only a small fraction (~ 1%) of the bolometric luminosity. 
In contrast, the regions emitting the UV and EUV com- 
ponents are located close to the BH, typically at distances 
of a few 10 Rq. Modeling the disc spectrum in the UV 
and EUV ranges is a difficult task as radiation is subject 
to several processes such as Comptonization. The emer- 
gent spectrum does depend strongly on the details of the 
disc vertical stratification which is not well known. On 
the contrary, in the optical range, one can show that the 



local spectrum is roughly a blackbody at the effective tem- 
perature corresponding to the conversion of gravitational 
energy into radiation (Collin 2000). 

The bolometric luminosity Lboi is a robust prediction 
for a (steady state) standard disc: it is related only to the 
accretion rate M and to the radius of the last, dynami- 
cally stable orbit, and corresponds to an efficiency of mass- 
energy conversion r\ w 10% (the value adopted throughout 
this paper), depending on the spin of the black hole. Thus: 
L bo i = nMc 2 , (1) 
and in principle, accretion stops if Lboi > LEdd where 
LEdd — 1.5 x 10 45 M7 erg/s is the Eddington luminosity 
and Mj is the mass in units of 10 7 M Q . 

AGN spectra also display a hard X-ray spectrum which 
extends up to a few hundred keV. This hard component 
which contributes a significant fraction of the bolometric 
luminosity (typically 30%) cannot be produced by a bare 
standard disc. This is why one generally invokes the pres- 
ence of a hot gas close to the BH (see Fig. |[ left panel), 
either in the form of a corona above the disc, or, instead 
of the inner disc itself, a quasi-spherical medium. In this 
case, the luminosities emitted respectively in the X-ray 
and optical-UV ranges relative to the accretion luminos- 
ity depend on the radius where this hot medium physically 
connects to the outer, standard disc. 

The likely change in the morphology of the accretion 
flow in the inner region (thin disc — > thick disc/spherical 
medium) is predicted by theory. The standard model prob- 
ably does not work at small radii where the disc is radi- 
ation pressure supported because of thermal instabilities. 
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In comparison, the standard disc as a model for outer re- 
gions (R ~ a few 10 2 -Rg) seems more consistent: the disc 
is sustained by gas pressure, relatively cold, and much 
more stable. It is true that a thermal instability, possibly 
recurrent (Siemiginovska, Czerny & Kostyunin 1986), is 
expected where hydrogen recombines/ionizes but such an 
instability is predicted to keep the disc in a geometrically 
thin configuration. 



2.2. The status of the outer, self- gravitating disc 

At large radii (roughly R > I0 3 i?G! see below), the stan- 
dard disc solution must be modified to account for the disc 
vertical self-gravity which may exceed the vertical com- 
ponent of the central gravity. The occurrence of vertical 
self-gravity is inescapable since it varies roughly as R~ 3 ^ 5 , 
that is less steeply than central gravity which varies as 
R- 3 (Shlosman & Begelman 1998; Hure et al. 1994; Hure 
1998). When the hydrostatic equilibrium is no longer con- 
trolled by the central object but by the disc, the surface 
density continuously increases whilst the disc gets thin- 
ner and thinner. But a self-gravitating disc is expected to 
be gravitationally unstable as soon as some criterion is 
fulfilled (Goldreich & Lyndell-Bell 1965; see also Toomre 
1964), usually this is Q = < 1, where fi is the Keple- 
rian angular velocity and p is the local mass density. 

The structure of the disc in the gravitationally unsta- 
ble region is not known. It may become marginally stable 
again, if the viscosity self-adjusts in some fashion (Lynden- 
Bell & Pringle 1974). Gravitational instabilities increase 
random (radial and vertical) velocities and could lead to 
the disc fragmentation into clouds (Pacynski 1968; Shore 
& White 1982, Shlosman & Begelman 1989). Following 
this scenario (see Fig. |], right panel), it is often assumed 
that these clouds are moving with a high velocity, and 
that their chaotic motion and mutual collisions provide 
the support against vertical gravity and keep accretion at 
the required rate (Begelman & Krolik 1986, Kumar 1999). 
The clumpy disc could become geometrically thick again, 
and even, at very large distance (~ 1 pc), it could con- 
nect to the dusty molecular torus invoked in the Unified 
Scheme of Seyfert nuclei (Antonnucci & Miller 1985). It is 
also possible that fragments in the unstable disc collapse 
to form protostars. Because these protostars can trap a 
large amount of gas from the disc, they should rapidly 
evolve into massive stars, and then give rise to outflowing 
gas through winds and supernovae explosions (Collin & 
Zalm 1999). In both scenarios one would thus expect the 
presence of a dense medium above the accretion disc, with 
an azimuthal velocity close to the Keplerian one, which 
could be identified with the BLR. 



2.3. Size and velocities within the BLR 

A major tool for studying the structure of the BLR is 
reverberation mapping, through the study of correlated 
variations of the lines and continuum fluxes (see the re- 
view by Peterson 1993). It has been used to determine the 
size of the BLR in several tens of AGN, and to correlate 
the mass of the central black hole with the luminosity of 
the AGN, under the assumption that the BLR is gravita- 
tionally bound to the black hole. 

Though powerful for size determination, reverberation 
mapping does not provide the velocity field of the BLR, 
and can generally not distinguish between radial and ro- 
tational motions. In some objects however, there are in- 
dications that rotational motions of the medium emitting 
the H/3 line dominates over radial motions. On the other 
hand, the medium emitting high ionization lines (or HILs) 
is sometimes observed as an outflowing medium. This is 
thought to be the case in NGC 4051 (Peterson et al. 2000). 
However it is likely that the velocity Vblr does not differ 
strongly from a gravitationally bound, Keplerian or virial 
motion (Peterson & Wandel 1999). Then, the radius of the 
BLR is approximately given by 
Rblk GM c 2 900 

R G Vi LR X GM ~ Vi [) 
where Vblr = Vq x 10 9 cm/s. 

Because there is no consensus about the dynamics of 
the BLR, there has been little progress in understanding 
the physics of this region. The most commonly accepted 
picture, deduced from photoionization models, consists of 
an assembly of photoionized clouds that cover at least 10% 
of the source emitting the incident, primary continuum. 
The number of BLR clouds is probably very large because 
the volumic filling factor of the emitting medium is ex- 
tremely small. This large population ensures a large total 
emitting area as well as a high density and a small individ- 
ual size. The smoothness of the line profiles also requires 
a large number of clouds. If these clouds are not confined 
by any mechanism, they expand in a very short time scale 
(less than one year). If the BLR is stationary, new clouds 
must be permanently generated. Alternatively, if the the 
clouds are formed uniquely then some form of confinement 
mechanism must operate. The thermal confinement by a 
hot medium as proposed in the "two phase model" by 
Krolik, Tarter & Mc Kee (1981) implies a Thomson thick 
hot medium, which is difficult to reconcile with the very 
short variation time scale of the X-ray flux. 

A natural idea which does not address the confine- 
ment problem is that the BLR is made of the atmosphere 
of giant or "bloated stars" , or by winds from giant stars 
(Edwards 1981, Periston 1988, Scoville & Norman 1988, 
Alexander & Netzer 1994). This model however requires 
a very large number of stars to account both for the line 
luminosity and for the smoothness of the line profiles. An- 
other possibility that does not appeal to any special con- 
finement mechanism either is that the lines are emitted by 
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the accretion disc. Dumont & Collin-Souffrin (1990) have 
proposed that low ionization lines (or LILs), like Balmer 
lines, which require a high density emitting medium (Collin- 
Souffrin et al., 1986) are formed at the surface of the ac- 
cretion disc. On the other hand, Murray & Chiang (1995, 
1997) proposed that the broad emission lines are emitted 
by a wind released at the top of the disc where motions 
are still mainly rotational. In this model, the broad ab- 
sorption lines (BALs) are produced by the high-velocity 
component of the wind. In both cases, the BLR is domi- 
nated by rotation and its presence (or at least, that of the 
region emitting the LILs) is linked to the accretion disc. 
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3. Relation between the optical luminosity and 
the bolometric luminosity 

3.1. Global investigation 

Kaspi et al. (2000) combined their own reverberation map- 
ping study with other published data to determine the size 
of the BLR, the BH mass and the optical luminosity. As 
it is usually done (Peterson 1993, 2000), they derived the 
size of the BLR from the measurements of the time lags 
between the light curves relative to the Ha and H/3 fluxes 
and the light curve of the underlying continuum which is 
assumed to follow the ionizing continuum. From the size 
of the BLR, they computed the BH mass from the FWHM 
of H/3 and Ha, assuming that these lines are emitted by a 
gravitationally bound medium. There is some uncertainty 
in the mass determination (by a factor of a few) since 
there are two different, non-equivalent methods to deter- 
mine the FWHM of lines. The FWHM can be averaged on 
all spectra, giving a "mean" BH mass: this is the method 
preferred by Kaspi et al. (2000). In the method proposed 
by Peterson et al (1998), the FWHM is measured on the 
rms spectrum, giving a "rms" mass. 

To estimate the bolometric luminosity from the monochro- 
matic luminosity, Lboi °c ^L° bs (5100A) is generally as- 
sumed. For the AGN generic spectrum (Laor et al. 1997), 
one deduces: 
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Fig. 2. Z^ isc (5100A)/L° bs (5100A) ratio versus the mass 
(in units of 10 7 M Q ) for the Kaspi et al. (2000) sample: 
L^ 1SC (5100A) have been computed using Eqs. (0) and (||). 
We have distinguished NLS1 galaxies (squares) and other 
nuclei (circles). 



9 x z/L° bs (5100A). 



(3) 



For each object in Kaspi et al. sample, we have compared 
this optical-to-bolometric conversion with two non equiv- 
alent but both self- consistent methods of conversion based 
on the disc model: 

1. by computing the theoretical, monochromatic luminos- 
ity i/L^5 lsc (5100A) from the spectrum of an optically 
thick, steady state disc given the pair (M, M), assum- 
ing Eq.(H), and Eq.(Q) to estimate M. Figure || displays 
the ratio L^ sc (5100A)/L° bs (5100A) obtained with this 
method. It shows that, for all objects in the sample, the 
emission due to the disc can represent only a minor, 
and in most cases, a negligible fraction of the observed 
luminosity in the optical band. 



by finding the "right" accretion rate M (and then 
Lboi via Eq.(|l|)) such that the observed monochro- 
matic luminosity is only due to a standard disc (that 
is 2^ isc (5100A) = L° bs (5100A)). This method leads to 
extremely small ratios ^L° bs (5100A)/Lb o i, as shown 
by Fig. ^|. It implies a very large bolometric luminos- 
ity, and therefore a very large Eddington ratio rh — 
Lboi/^Edd- To illustrate this last point the correspond- 
ing Eddington ratio is plotted versus the central mass 
on Fig. ^. It shows that almost all the objects would 
be radiating at highly super- Eddington rates, proving 
that this second self-consistent optical-to-bolometric 
conversion is wrong. 

In other words, a relation of the form Eg. (^) should be 
roughly valid, meaning that the optical luminosity is not 
radiated by a standard disc. 

Note that we assume a face-on disc, thus overestimat- 
ing slightly L^ lsc . These conclusions are unchanged if the 
inner disc is removed for some reason, that is, if the inner 
radius is larger (in the limit of a few 10 Rq) than a few 
Rq as considered here, since the optical emission is almost 
entirely produced in outer regions located much further 
away from the center. Clearly the situation becomes even 
more extreme if the region emitting the optical band is 
not present in the disc, that is, if the disc is truncated for 
some reason, for instance due to self-gravity. 

There are at least four possible explanations for this 
result: 

1. the optical luminosity is not due to the disc, as quoted. 
This is plausible since the standard disc spectrum is 
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Fig. 3. L(5100A) o bs/£boi ratio versus the mass (in units 
of 10 7 M Q ) for the Kaspi et al. (2000) sample: L hol 
has been computed assuming that the optical luminosity 
L(5100A) o bs is entirely due to the disc. We have distin- 
guished NLS1 galaxies (squares) and other nuclei (circles). 
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Fig. 4. iboi/iEdd ratio versus the mass (in units of 
1O 7 M0), for the Kaspi et al. (2000) sample, assuming that 
the optical monochromatic luminosity L obs (5100A) is en- 
tirely due to the disc. We have distinguished NLS1 galax- 
ies (squares) and other nuclei (circles). Note the point cor- 
responding to the lowest mass is NGC 4051. 



well known to be flatter than the spectrum of Seyfert 
1 nuclei in the optical range (Koratkar & Blaes 1999, 
Hubeny et al. 2000, Collin 2000). This is also compat- 
ible with a recent study of the NLS1 1 RE J1034+396 
by Puchnarewicz et al. (2000) who conclude that an 
additional underlying power law connecting the IR to 
the X-rays is really necessary to fit the optical/UV 
continuum. They suggest that this extra component 
might be a non thermal BLac-type component. This is 
indeed possible if for instance there is a strong mag- 
netic field anchored into the disc. But it could also be 
due to a thermal process like the emission of dense 
clouds surrounding the X-ray source and reprocessing 
a large fraction of the X-ray luminosity, as proposed by 
Celotti, Fabian and Rees (1992), and by Collin-Souffrin 
et al. (1996). This cloud system should be distributed 
quasi-spherically in order to reprocess a large fraction 
of the X-ray luminosity (in the case of a disc, the repro- 
cessed fraction decreases as like the accretion lu- 
minosity, and is necessarily small in the optical range), 
the standard disc model is not valid. Although the ac- 
cretion luminosity does not depend on viscosity, the 
fact that the disc is effectively optically thick and ra- 
diates like a black body is not independent of the vis- 
cosity prescription. Note, as a proof, that the opti- 
cal thickness of a standard disc is a function of the 
a-viscosity parameter. Either a completely different 
mechanism of angular momentum evacuation in the 
disc is at work (e.g. density waves), or the formalism 
for turbulent viscosity must be reconsidered, or both, 
the accretion rate is not uniform. This possibility in- 
cludes the presence of a non steady disc and the possi- 
ble existence of some steady inflow involving an accre- 
tion rate decreasing with decreasing radius. Concern- 
ing the first point, there are a few well known sources 
of instability that can perturb the equilibrium of a disc 
but this would mean that all objects are observed in 
a "high optical state". This is statistically somewhat 
doubtful. Regarding the second point, this would mean 
that a huge amount of gas is accreted at large distances 
and then converted into an outflow at smaller radii. 
This is also unrealistic because it would imply massive 
super-Eddington outflows, which are not observed, 
accretion is advection dominated in the inner regions. 
Such a solution is still not demonstrated on theoretical 
grounds. Stable Advection Dominated Accretion Flows 
(ADAFs) (Narayan & Yi 1995, and many other works) 
take place for relatively low accretion rates, well below 
the Eddington limit. At the other extreme, there are 
advective solutions for Eddington rates ("slim discs" 
for instance, cf. Abramowicz et al. 1988) or for super- 
Eddington rates (thick discs), but it is difficult to con- 
ceive a durable, super-Eddington accretion rate, as it 
would imply a dramatic growth of the hole over cosmic 
time scales. In conclusion, even if the disc becomes ge- 
ometrically thick in the inner regions, advection is not 
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iboi/i obs (5100A) is larger in NLSls than in BLSls. How- 
ever this is impossible, as it would imply super-Eddington 
rates. Indeed Fig. |], which displays -Rblr / -Rg as a func- 
tion of the Eddington ratio, for L bo i/£ obs (5100A) = 9, 
shows that with the conservative assumption of Eq. [|, 
NLSls are already radiating at Eddington rates. NLSls 
occupy a privileged position on this figure, not only be- 
cause of their small FWHM (corresponding by definition 
to large -Rblr/-Rg values), but also because they display 
large Eddington ratios. Actually this was strongly sus- 
pected on the basis of their spectral and variability prop- 
erties in the X-ray range, similarly to galactic BH candi- 
dates (Pounds, Done & Osborne 1995). 

We conclude that NLSls have not only larger Edding- 
ton ratios, but also probably a larger "non disc" fraction 
of optical emission, compared to BLSls. 



Fig. 5. The points give -Rblr in lt-days versus the 
monochromatic luminosity at 5100 A in units of 10 44 erg 
s _1 , for the Kaspi et al. (2000) sample. The big open 
squares show NLSls. The curves give -Rcrit from the 2D 
simulations. The accretion rates of the theoretical curves 
have been converted into -L5100 assuming that the bolo- 
metric luminosity is equal to 9xL5100. The curves are 
labeled with the black hole mass in log(M Q ), and the Ed- 
dington ratio m is indicated on the curves. 

likely to be important in the Seyfert galaxies and mod- 
erately luminous quasars considered here, and the gas 
should flow all the way down to the last stable orbit 
and radiate accordingly. 



3. 2. The case of Narrow Line Seyfert 1 

Figure [5] displays the size of the BLR, -Rblr, versus vL^ lsc (5l 
for the Kaspi et al. (2000) sample. In this figure and in 
Figs. [| [| and[| NLS1 have been isolated^. We notice that 
they do not occupy any privileged region on Fig. ||, con- 
trary to Figs. [| H and ||. On Fig. Q they are concentrated in 
the top-left region corresponding to small masses and high 
bolometric luminosities. Further, NLSls have the smallest 
^ isc (5100l)/L° bs (5100A) ratios compared to other ob- 
jects in the same range of BH mass, as shown by Fig. g. 
How to explain these properties? 

We know that NLS1 galaxies have a strong soft X-ray 
excess compared to Broad Line Seyfert 1 galaxies (BLSls) 
(Boiler, Brandt & Fink 1996 for a review; see also the 
proceedings edited by Boiler et al. 2000). This would sug- 
gest that the BBB is shifted towards the EUV, and that 
the ratio L disc (5100l)/L disc (EUV) is smaller in NLSls 
than in BLSls. Meanwhile, this would mean that the ratio 



4. What determines the size of the BLR ? 

4-1. Model for the outer disc 

Based mainly on distance arguments, we have suggested 
before that the BLR is possibly related to the suppres- 
sion of the disc by gravitational instabilities. We have 
investigated this idea more quantitatively by a series of 
bi-dimensional simulations of steady state, Keplerian ac- 
cretion a-disc (Hure 2000). More precisely, the model in- 
cludes convection in the framework of the Mixing Length 
Theory, turbulent pressure that substantially thickens the 
disc and self-gravity within the infinite disc approxima- 
tion. Opacities and equation of state are realistic and cor- 
responds to a gas at LTE with cosmic abundances. Ex- 
ternal irradiation is not taken into account, for several 
reasons. First the discs considered here are optically very 
thick and moreover the irradiation flux is likely to decrease 
le gravitational flux, having thus a relatively small in- 
uence on the disc structure (Collin & Hure 1999). Sec- 
ond, as we shall see below, outer regions are in the shadow 
of the central regions due to self-gravity which pinches 
the disc vertically. Figure ^ shows an example of internal 
structure of a disc computed for Al = 10 s M Q , M = 10 _1 
M Q /yr and a = 0.1. We see that the disc becomes self- 
gravitating (in the sense AirGpR? > M) from ~ 1200 Rq 
which results in a significant decrease of the disc thickness 
and rise of the density. Note however that the disc flaring 
vanishes at much lower radius ~ 900 Rg. 



1 We define NLS1 as Seyfert nuclei with Vblr = 
(GM/ -Rblr) 1/2 < 2000 km/s, instead of using FWHMs, to 
avoid having to take into account a correction for Ha and H/3 



4-2. Location of the gravitationally unstable region 

In order to determine the radius -R cr it where the disc be- 
comes gravitationally unstable, we have computed a series 
of disc models by varying M and M. There are many pos- 
sible sources of uncertainties on the value of this quantity, 
including a lack of information regarding some physical 
processes like turbulent viscosity through the cv-prescription 
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Fig. 6. Density map (top) and temperature map (bottom) 
for a steady state Keplerian a-disc, surrounding a 10 8 M 
black hole, versus R and z both normalized to Rq. The 
accretion rate is M = 1CP 1 M Q /yr (i.e. m ~ 0.03), and 
a = 0.1. 

and its z-dependency, convective transport, self-gravity, 
environment effects. In addition, the threshold value of 
the Q-parameter which defines the instability criterion is 
uncertain in this context (e.g., it is not clear whether this 
criterion applies to the midplane quantities). We have at- 
tempted to estimate an absolute error on i? cr it by consider- 
ing 0.01 <a<landi<Q<4 (Toomre, 1964; Goldreich 
& Lynden-Bell 1965). Figure |^ shows the location of the 
unstable disc region (Rent ~ 1500 — 2000 Rq)- Our results 
are displayed in Figure which gives Rent versus the ac- 
cretion rate for four mass decades, obtained for a = 0.1 
and Q = 1 (mean values). We see that i? C rit is mainly de- 
termined by the central mass and weakly sensitive to the 
accretion rate M. The a-parameter and threshold value 
for Q have rather minor effect for moderate and high Ed- 
dington luminosities. The magnitude of the error on i? cr it 
is given by the gap between curves corresponding to the 
same central mass but to different values of a and thresh- 
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Fig. 7. i? C rit for Q = 1 as a function of the accretion rate 
expressed in Mq/yt. The curves are labeled with the black 
hole mass in log(Mo). a — 0.1 throughout. A few other 
cases are also displayed (for a = 0.01 and a = 1, and for 
= 4). 

old values for Q. It is comforting to see that the results 
are quite reliable, as the uncertainties do not produce an 
"error" larger than a factor say, ~ 3 on i? cr i t . 

Figure ^ displays R cr it/RG versus the Eddington ratio 
for the same models. A crude fit in the domain 0.01 < 
rh < 1 gives 

^ w2 x 10 4 M 7 -°- 46 (4) 



4-3. Comparing the BLR and the disc sizes 

We have plotted on Fig. ||i?crit in light-days versus L(5100A), 
using Eq. (^) to convert the disc accretion rate into op- 
tical luminosity (otherwise, accretion rates are unrealis- 
tically high), together with the data of the Kaspi et al. 
sample. Figure g displays Rblr/Rgi for the Kaspi et al. 
sample, still assuming Eq. (||). We see clearly from these 
two figures that, in absolute, the BLR size is correlated 
with the occurrence of the disc gravitational instability. 
So, if we identify the edge of the gravitationally unstable 
disc with the inner edge BLR, that is i? C rit < Rblr, we 
conclude that active nuclei hosting low mass black holes 
are predicted to have broad lines narrower than others. 

Fig. H shows that, though the observed masses do not 
match exactly the theoretical expectations, there is a clear 
trend that objects with small masses (10 6 < M < 10 7 M©) 
have a large i?BLR/-Rcrit ratio and that objects with large 
masses (M > 10 s M Q ) have a small -RBLR./-Rcrit ratio, in 
agreement with our theoretical predictions. So, at least to 
this "zero order approximation" , the gravitationally un- 
stable region of the disc has kinematic and size properties 
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Fig. 8. The points give -Rblr/^g versus L ho i/L Edd , 
binned according to the BH mass, and assuming Lboi = 
9 x L(5100A), for the Kaspi et al. (2000) sample, circles: 
M < 10 7 M Q ; triangles: 10 7 < M < 10 s M ; squares: 
M > 10 8 M Q . The big open squares are NLSls. The curves 
give i? cr it/i?G from the 2D simulations. They are labeled 
with the black hole mass in log(MQ). 

compatible with the BLR and could therefore be the ul- 
timate source of the BLR clouds. It is however difficult 
to draw more quantitative conclusions from this investi- 
gation. 



4-4- Influence of the ionization parameter of the BLR 

One can notice that -R cr it depends mainly of the BH mass 
and little on the luminosity, while a correlation appears 
on Fig. | between .Rblr and the J L obs (5100l). Kaspi et al. 
(2000) represent this correlation as -Rblr oc L(5100l) a7 °. 
A smaller power law index (closer to 0.5) was found in 
previous works (e.g. Wandel et al. 1999). Anyway it is 
not clear how to translate such a relation into a relation 
between -Rblr and Lboi, because of an important uncer- 
tainty arising from the optical-to-bolometric conversion. 
If the optical to bolometric luminosity ratio is constant, it 
results in -Rblr oc L^ 01 . 

In a recent paper, Nicastro (2000) proposed that the 
BLR clouds are released by the accretion disc in the region 
where a vertically outflowing corona exists, according to a 
model proposed by Witt, Czerny & Zicky (1997). This is 
indeed an interesting model, since in this case the main pa- 
rameter governing the size of the BLR would be rh and not 
M . However one would expect that the remote regions giv- 
ing rise to the BLR, which are strongly pressurized by the 
corona, are also gravitationally unstable. Moreover there 



is also another important parameter governing the size of 
the BLR, that is its ionization parameter 



where L- lorl is the ionizing luminosity (~ £boi) and n is 
the number density in the clouds. Actually, any model 
must reproduce not only the kinematics and dimension of 
the BLR but also the ionization parameter (linked to the 
ionization state) consistent with the observed line ratios 
which are almost independent on the luminosity over more 
than four decades (except for the Baldwin effect which 
might indicate a decrease of the ionization parameter with 
increasing luminosity, but only for high ionization lines). 
Typically, £ ps 1 and detailed studies of the line ratios show 
that n pa 10 10 cm' 3 . If the density is about a constant, 
the constancy of £ translates into a relation -Rblr oc ^bol' 
in close agreement with the observations. Another way to 
express this relation is that the flux irradiating the clouds 
is constant. This relation has been used in the "ioniza- 
tion method" aiming to determine the mass of the black 
holes, and it gives results compatible with the reverbera- 
tion mapping method (Wandel, Peterson & Malkan 1999), 
but more uncertain. 

Thus a possible interpretation of the observed correla- 
tion is that among clouds released by the gravitationally 
unstable disc, only those located at the right distance (i.e. 
corresponding to the right ionization parameter or ion- 
izing flux) are observed as BLR clouds: the others could 
contribute to the Warm Absorber. This could explain why 
-R cr it constitutes a "lower envelope" to -Rblr for a given 
BH mass. Also since a constant ionizing flux corresponds 
to -Rblr/-Rg oc y/m/M, we see that the objects display- 
ing the largest -Rblr/^Rg ratios, hence the smaller line 
widths, are those which have the largest accretion rates 
in Eddington unit and the smallest masses. Such a situa- 
tion is proposed to explain the properties of NLSls. Note 
that in this picture, the more intense soft X-ray excess 
displayed by NLSls compared to BLSls can also play a 
role, since such a spectrum has a stronger ionizing power 
for the same ionization parameter (e.g. Wandel 1997). 

5. Conclusion 

We have discussed the relations between the optical lu- 
minosity, the BH mass, and the size of the BLR, in the 
framework of the standard accretion disc model. 

We have shown first that the optical luminosity can- 
not be accounted for, meaning that the standard accretion 
disc picture (stationary, geometrically thin and optically 
thick) does not hold, at least in the region emitting the 
optical band, i.e. at ~ 100-Rq. Either the major fraction 
of the optical luminosity is not due to disc emission, or 
the disc is not "standard" : it is unstable, or the accretion 
rate depends on the radius, owing to strong outflows or 
to advection dominated accretion. In these last cases, the 



Suzy Collin and Jean-Marc Hure: The Size-Mass-Luminosity Relations in AGN. 



9 



implied mass rate would have to be strongly super-critical. 
We would thus favor the first explanation. We have also 
shown that NLSls are extreme in this context: either they 
have a larger fraction of non disc emission in the optical 
range and are radiating close to their Eddington lumi- 
nosity, or they have strongly super-Eddington luminosity, 
which seems implausible. So again in this case it seems 
that the best explanation would be that a large fraction 
of the optical emission (larger even than for BLSls) is not 
produced by the disc. 

In the second part, we have studied the relation be- 
tween the disc and the BLR, and we have shown that there 
is a good agreement between the size of the BLR and the 
critical radius at which an a-disc becomes gravitationally 
unstable. This suggests that the BLR is produced above 
the gravitationally unstable part of the disc. However the 
ionization parameter should also play a role in explaining 
the correlation observed between the size and the luminos- 
ity. If BLR clouds are seen only for a small range of ion- 
ization parameter or of ionizing flux, the size of the BLR 
expressed in Rq should increase with the accretion rate 
expressed in Eddington unit and decrease with the BH 
mass, according to the observed correlation. This could 
explain the small widths of NLSls. 
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